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Abstract. We simulate macroscopic shear experiments in active nematics and compare them with mi- 
crorheology simulations where a spherical probe particle is dragged through an active fluid. In both cases 
we define an effective viscosity: in the case of bulk shear simulations this is the ratio between shear stress 
and shear rate, whereas in the microrheology case it involves the ratio between the friction coefficient and 
the particle size. We show that this effective viscosity, rather than being solely a property of the active 
fluid, is affected by the way chosen to measure it, and strongly depends on details such as the anchoring 
conditions at the probe surface and on both the system size and the size of the probe particle. 
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1 Introduction 

Active fluids are suspensions of particles that absorb en- 
ergy from their surroundings (or from an internal reser- 
voir) in order to do work. They represent a novel type 
of out-of-equilibrium system and they are found predom- 
inantly, although not solely, in biological contexts [1 . Ac- 
tive particles exert forces on the fluid in which they are 
embedded and, to lowest order, they can be modelled as 
a force dipole. Extensile active particles exert forces that 
are directed outward, along their main axis, while con- 
tractile particles exert forces directed inward. Examples 
of the former are suspensions of bacteria like E. coli, while 
dispersions of algae, such as Chlamidomonas^ and the ac- 
tomyosin solution which constitutes the cytoskeleton of 
eukariotic cells, are contractile fluids. 

Due to the continuous stirring exerted by their con- 
stituent particles, active fluids show a wide range of novel 
properties, the most striking of which is probably the ap- 
pearence of a spontaneous flow phase: for high enough 
dipolar forces or active particle density, these fluids can 
flow spontaneously in steady state, in the absence of any 
external force [lEHEllS] . 

The bulk rheological properties of active fluids are also 
highly non-standard. In particular, theory predicted |7l 
IH], and both simulations [9 and experiments [10,TT1 con- 
flrmed, that active nematics increase their viscosity with 
activity if contractile and decrease it if extensile. On the 
other hand, microrheology experiments, where the motion 
of micron-sized probe particles is studied [12], allow one 
to investigate locals rather than global, properties which 
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are important in biophysical contexts. (For instance, cells 
are often subjected to localised cues rather than to bulk 
external body forces These experiments have also 

shown strong violations of the fluctuation-dissipation the- 
orem [HITS], further proving that the physics of active 
fluids has little in common with that of passive, near- 
equilibrium systems. 

Here we present a systematic numerical study of the 
rheological properties of active fluids, which we investigate 
both globally, through shear experiments (in a quasi-one 
dimensional geometry, as used in [9 ), and locally, by ap- 
plying a constant force on a spherical colloid embedded in 
the fluid. The resulting comparison is new, and important, 
as it allows us to probe to what extent an effective viscos- 
ity, ^eff, exists for active fluids: for this concept to be a use- 
ful one, T^eff should not depend much on the way in which 
the behaviour of the active fluid is probed. It hence should 
have the same or a similar value in both bulk rheology and 
microrheology experiments. For bulk rheology simulations 
the effective viscosity is defined as the ratio between the 
steady state shear stress and the applied shear rate. For 
microrheology, if we measure the steady state velocity v 
reached by a spherical particle subjected to a pulling force 
F, the effective viscosity can be defined as r^eff = ^o/^Tri?, 
with ^0 = F/v the friction coefficient. 

Our bulk rheology simulations show that in sheared 
nematics ?7eff differs for different anchoring conditions at 
the container walls, as expected in view of the anisotropic 
nature of the fiuid. However, the effective viscosity also 
depends on system size, in a way which it is possible to 
predict very accurately from an analysis of the underly- 
ing equations of motion. With our microrheology simula- 
tions, instead, we find that, when a spherical particle is 
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dragged through an active fluid, the drag force does not 
depend linearly on the particle radius, i?, in strong viola- 
tion of Stokes' law. We first present our results in the case 
of planar anchoring at the probe surface and we explain 
this effect in terms of the active stress distribution and 
the director-field deformation determined by the anchor- 
ing condition on the particle surface. To test this reasoning 
we consider also the case of normal anchoring and the ex- 
treme case in which no anchoring condition is imposed, so 
that the far field order is kept also in the proximity of the 
particle. 

While we performed most microrheology simulations 
with periodic boundary conditions, we also checked the 
role of boundary walls. To this end, we repeated the mi- 
crorheology experiments relative to planar anchoring at 
the probe surface, when the sample was bounded by walls 
along one direction, and found good qualitative agreement 
with the case of periodic boundary conditions. Finally, in 
order to connect more easily to future experiments, we 
modelled also the case in which the particle disrupts the 
active fluid ordering, by imposing zero order parameter on 
its surface. We did this because the particles used in mi- 
crorheology experiments are usually micron-sized, so that 
the size of the microrheological probe may be in practice 
similar to that of the active particles (e.g. suspensions of 
bacteria, such as E.coli). In this case one may envisage 
that, rather than creating homogeneous/normal anchor- 
ing, the probe particle may simply locally create a more 
disordered region. 

Our paper is structured as follows. In Section [2] we 
present our model. In SectionjS we discuss our results, first 
for bulk rheology (Section 
on microrheology (Section 



3. ID and then for simulations 
3.2). In Section |4] we present 



further discussion on our results and draw conclusions. 
Some elements of the microrheology work were reported 
briefly in [16], which mainly focuses however on a regime 
of high activity /large probe size that is not addressed here. 



2 Model and methods 

In the continuum limit, active fluids can be described 
within a phenomenological hydrodynamic model, which 
is a generalisation of the equations of motion describing 
liquid crystal hydrodynamics. The additional terms which 
we introduce account for the effect of forces exerted by the 
active particles on the fluid. Liquid crystals are often de- 
scribed in terms of the director n, a head-tail symmetric 
vector describing the local average orientation of the rod- 
like particles which typically constitute them. A more gen- 
eral characterisation can be achieved by defining a tenso- 
rial order parameter Q, which can account for the average 
magnitude of orientational order, on top of its direction. 
The magnitude of order, and the director field, n, are 
then associated to respectively the largest eigenvalue of Q 
and its associated eigenvector. 

In this work, we start from a generic Landau - de 
Gennes free energy density T — Th ^ J^ei + which 
describes the equilibrium physics of liquid crystals and 



consists of a bulk, of an elastic, and of a surface contri- 
bution. The bulk contribution, in the nematic case, takes 
the following form [17 : 
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(1) 

where Aq is a constant, with dimension of free energy den- 
sity, while the parameter 7 controls the magnitude of order 
- this is related physically to either the temperature (in 
thermotropic liquid crystals), or the concentration (in ly- 
otropic liquid crystals). The transition from the isotropic 
to the nematic phase takes place at 7 = 2.7; here we work 
at 7 = 3.0, which leads to an ordered phase. The elastic 
free energy density is: 



K 



(2) 



where K is an elastic constant (we fixed it to i^^ = 0.05 
in simulation units (s.u.) in what follows, see below for a 
mapping between s.u. and physical units). The one elas- 
tic constant approximation which we choose here does 
not modify the physics we are interested in qualitatively 
(for an example of the effect of using different splay and 
bend elastic constant in active fluids, without colloids, see 
e.g. [18 ). Note that in Eqs. ^ and ([2|, and below, Greek 
indices are used to denote Cartesian coordinates, and that 
summation over repeated indices is implied. 

The surface free-energy density term that describes the 
anchoring of nematogens to a solid surface (either a wall 
or a colloidal surface) is: 



(3) 



where Q^^ = 5*0 (n^n^ — Sa/s/^), So quantifies the pre- 
ferred degree of order, and describes the director ori- 
entation on the surface [19]. Different types of anchor- 
ing are implemented in our code by fixing accord- 
ingly (see ^20f j21 ^ . e.g. tangential anchoring is enforced 
by projecting the order parameter on the local colloidal 
surface [22^), whereas W is the parameter controlling the 
anchoring strength. The boundary conditions for the or- 
der parameter on the colloidal surface are given by [23l 
HU: 

where is the local normal to the colloid surface. 

The hydrodynamic equation for the evolution of the 
order parameter is: 



(S, + n- V)Q-S(Vn,Q) = rH, 



(5) 



where is a collective rotational diffusion coefficient re- 
lated to the rotational viscosity of the liquid crystal [2[|3l 
4 . The second term of eq. ^ describes the advection of 
the order parameter; while S(Vu, Q) describes the possi- 
bility for Q to be rotated and stretched by the fluid and 
takes the form [25l: 



S(Vn,Q) = (eD + a;)(Q + l/3) + 

+(Q + l/3)(^D -uj)- 2^(Q + l/3)Tr(QVn). (6) 
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In Eqj6] Tr denotes the tensorial trace, while D = (Vn + 
Vu^)/2 and uj = {yu — Vu^)/2 are respectively the sym- 
metric and the antisymmetric parts of the velocity gradi- 
ent tensor. The right-hand side of eq. (5| describes the 
relaxation of the order parameter towards an equilibrium 
configuration. The tensor H {molecular field) represents 
the symmetric and traceless part of the functional deriva- 
tive of the free energy F = J d^xJ^{x) with respect to the 
order parameter Qais- 
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The evolution of fluid momentum for an incompressible 
fluid {daU^s = 0) is described by a modified Navier-Stokes 
equation where new stress terms are introduced: 

p{dt^U(3d(3)Ua = Vd(3{d(3Ua^daUf3)^df3{n^^^n^^). (8) 

Here 77 is the bare viscosity of the Newtonian fluid (set 
to 0.6 in S.U.), while 77^^ is the stress tensor of a passive 
liquid crystal, whose detailed expression is [25] : 

daQju ~^ QajHji^ Hq^^Q^^ ^ (9) 

where Pq is the isotropic pressure. The parameter ^ deter- 
mines whether the liquid crystal is flow tumbling (which 
happens for ^ < 0.6) or flow aligning (^ > 0.6). In the 
absence of activity but with flow present, particles in flow 
aligning liquid crystals assume a fixed angle with respect 
to the flow direction (Leslie angle), while in flow tumbling 
they continuously change their orientation, chaotically. To 
avoid further complexity arising from this, we fixed ^ = 0.7 

(S.U.). 

Hydrodynamic equations for active fluids are obtained 
through an additional stress term in eq. |8| which takes 
into account the effect of active forces, and is given by: 



(10) 



where the activity coefficient, is a constant related to 
the size of active particles, to their density and to the 
intensity of the forces they exert on the ambient fluid. 
The parameter C, is negative for contractile and positive 
for extensile liquid crystals [T. 

We solve these equations through a hybrid Lattice- 
Boltzmann method (described in [201121] ) , where eq. [s] 
is solved through a finite-difference algorithm, while the 
Lattice-Boltzmann method is used for momentum trans- 
port (Navier-Stokes equation and continuity). Within this 
framework, a spherical probe colloid is simulated through 
the method of bounce-back on links, which allows us to 
implement the no-slip boundary condition at the colloid 



surface [26] . Different stress terms appearing in eq. [Sj have 
to be integrated over the colloid surface to determine the 
force and the torque acting on the colloid due to the sur- 
rounding fluid [20l[2ll[23]. 

Above and in what follows, parameters and results are 
presented in simulation units. To convert them into phys- 
ical ones, relative to a contractile actomyosin solution, we 
can consider values holding for typical intracellular actin 
gels ([27 ): the elastic constant K is 1.25 nN, and the 
rotational viscosity is 10 kPa/s. In this way (see [SjElH] 
for similar mappings) one simulation unit for forces cor- 
responds to 25nN, the lattice unit for length, Z\x, is 0.5 
/im, while the time unit corresponds to 10 ms [^. Note 
that the mapping to a different active fluid, such as a bac- 
terial suspension, would be significantly different [2 ,3, 4 . 
Finally, the anchoring strength W ^ when different from 
zero, is set to Kj Ax. 

Note that the above equations assume nematic order, 
in which individual active particles, which generally are 
polar, are equally likely to point in any direction along 
n. For polar phases different equations are needed (see 
e.g. [30]). Eq. (10) however holds at first order for both 
polar and apolar phases (see [r,7]), and we therefore ex- 
pect rheological properties not to strongly depend on po- 
larity. We will thus restrict ourselves to the apolar case in 
this work. 

We also point out that, as we employ a continuous 
model for active liquid crystals, our results on microrhe- 
ology are reliable only as far as we address colloids with 
a size much larger than the size of the active nemato- 
gens, e.g. the length of a bacterium, or the mesh size of 
an actomyosin solution. Typical probes in microrheology 
are in the 1-10 /im range. Therefore in order to describe 
bacterial or algal microrheology within the continuum ac- 
tive gel theory, one may need to use rather large probes. 
However this restriction is lifted in the case of contractile 
actomyosin gels. 



3 Results 

3.1 Macrorheology 

Bulk rheological properties of active fluids have been in- 
vestigated, theoretically [7l|8], with simulations [9]l31 . and 
through experiments [TOl[TT]. 

Given the importance of boundary conditions in de- 
termining the behaviour of liquid crystalline and of active 
fluids (as found e.g. in [UHllQ]), we address here the ques- 
tion of the role played by the anchoring condition on the 
walls and by the system size. With this aim in mind, we 
consider two different cases: planar anchoring, where ac- 
tive nematogens are oriented along the walls, and normal 
anchoring. 

As anticipated in the Introduction, an effective viscos- 
ity ^eff can be defined as the ratio between the measured 
stress and the imposed average shear rate, and we measure 
it for different values of C,. We highlight that here we are in- 
terested in the bulk microrheology of active fluids in their 
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Fig. 1. Dependence of ?7eff on system size L in the case of contractile active nematics with planar (a) and normal (b) anchoring, 
and an extensile active fluid with planar (c) and normal (d) anchoring at the walls. Legend on the left refers to both graphs in 
the same line. The inset in (d) shows the behaviour for ^ = 0.0003 more in detail. Points here refer to numerical results, while 
solid lines are theoretical predictions obtained through Eq. |18| 



quiescent (passive) phase, i.e. when they are not sponta- 
neously flowing. We focus here for simplicity on the case 
of a quasi-one-dimensional system, where there is transla- 
tional invariance in the vorticity and flow direction, and 
all quantities only vary along the flow gradient direction, 
which is z in our notation. This case is already sufficient 
to describe qualitatively the basic features of the spon- 
taneous flow transition in active fluids f2]|3l[4l|9]. Sponta- 
neous flow can be avoided by applying a small enough 
stress. 

Fig. [l] shows that 77eff in the nematic phase in our ge- 
ometry behaves consistently with previous theoretical ex- 
pectations [7,8 : it is higher with respect to the passive 
isotropic viscosity in the case of contractile active fluids 
and smaller in the extensile case. Moreover, by analysing 
its dependence on the system size L, which we define as the 
distance between the two shearing planes, we find that ryefr 
increases with L in contractile active fluids and decreases 
with L in extensile ones. As can be seen from Fig. [T] the 
anchoring of the active particle orientation at the walls 
plays quantitatively an important role: for instance the 



difference with the passive viscosity is enhanced for nor- 
mal anchoring. 

An analytical prediction for the apparent viscosity r^eff 
at low shear rat es is obtained by linearising the equations 
of motion, Eqs. |5|8[ Assuming that the flow is in the y- 
direction, we obtain 



rKSQ'^,{z) + ^u'y{z)=0, 



0, 



(11) 
(12) 



where SQyz{z) and Uy{z) are the deviations of the order- 
parameter tensor and the velocity from their rest values 
due to the applied shear; prime denotes the derivative with 
respect to z. Here, 



a = ^ — 1 X = 
for the planar anchoring, and 

a = ^ X = 



2e + - 3) 



2C + 5(^ + 3) 



(13) 



(14) 
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for the normal anchoring. The linear reponse of the system 
to the applied shear, satisfying the boundary conditions 
Uy{0) = -vo, Uy{L) = vq, and ^^^^^(O) = SQyz{L) = 0, is 
given by 



Uy{z) = VO 



SQyz{z) 



si nil Xz 

sinh ^ ' 



avo 



2rK\ tanh 



AL 



cosh \z 
cosh 



where 



A 



2r]rK + Kx'^ ' 



(15) 



(16) 



(17) 



The apparent viscosity is calculated by dividing the total 
shear stress by the shear rate 2vq / L which yields 



^eff 



LA 



2 tanh 



XL 



X 

2r 



(18) 



In Fig. [T] we compare this analytical prediction (solid 
lines) with our lattice Boltzmann data (symbols). The re- 
sults show a very good agreement between the viscosity 
in the linear regime obtained in our simulations with the 
predictions in Eq. 18, except for the case of an extensile 
fluid with L = 100 and ( = 0.0003, where spontaneous 
flow sets up in the numerical solution at 7 = 2 • 10~^, 
which introduces non-linear effects invalidating our lin- 
earised analytical treatment. 

A key feature of our results is that the effective vis- 
cosity for the active (but not passive) system strongly de- 
pends on system size, so that this quantity would be ill 
deflned in a system where the walls are inflnitely far apart, 
clearly at odds with what is expected of a "well-behaved" 
macroscopic viscosity of a fluid sample. In other words, 
the viscosity of an active fluid cannot be operationally de- 
flned as the ratio between shear stress and shear rate, as 
routinely done in experiments probing the bulk rheology 
of complex fluids. Notably, as Fig. [T] shows, this is despite 
the existence of a linear regime in the flow curve, which 
does appear and can be fltted with our analytical formu- 
las. 



3.2 Microrheology 

We now turn to the microrheology of active fluids. To 
study this, we performed simulations in which a spheri- 
cal particle is dragged by a constant force through an ac- 
tive nematic, considering the cases in which the dragging 
direction is either parallel or orthogonal to the average 
direction of the nematogens in the bulk, which we will re- 
fer to as the far fleld director. These studies simulate the 
simplest class of microrheological experiments: for passive 
fluids in the linear regime, where the drag force on the 
particle depends linearly on its velocity, one would expect 
to see Stokes' law at work. This prescribes that the drag 
coefflcient should depend linearly on particle size and on 
the fluid viscosity. Stokes' law is indeed typically used in 



experiments to infer mean properties of the fluid, in par- 
ticular its viscosity. 

Here we flrst monitor the colloid velocity in steady 
state, V, as a function of the magnitude of the applied 
force, F. From Fig. [2] it can be seen that these two quan- 
tities are linearly related, thereby allowing us to deflne 
the drag on the particle as ^0 = F/v. We then deflne an 
effective viscosity via Stokes' law as 77eff = ^o/^Tri?. 

We performed most of our simulations with periodic 
boundary conditions, for which corrections taking into ac- 
count hydrodynamic interactions between different images 
of the same particles are needed. These are calculated in 
the case of a simple fluid in [32 , while to our knowledge 
no similar calculations were done for liquid crystals. How- 
ever, we found that [32 seems to work quite well also in 
this case, since corrected data for passive nematics do not 
show signiflcant dependence of r^eff on R. We then applied 
the same corrections in the case of active fluids. 

Just as we considered different system sizes, L, in the 
macrorheology section, here we study the behaviour for 
different particle sizes R. 

We chose values of ( high enough for activity to af- 
fect the system behaviour, but small enough not to see 
any spontaneous bulk flow effects. In particular, we con- 
sidered ( = —0.001 and —0.002 for contractile nematics 
and C = 0.0001 and 0.0002 for extensile fluids [33], in the 
case of planar anchoring. In all other cases we chose a sin- 
gle value for C (C = —0.001 for contractile and ( = 0.0001 
for extensile nematics), as we wish to focus on the effect 
of anchoring on r^eff- 

We flrst analyse the case when the anchoring is planar, 
i.e. n lies in any direction on the plane which is locally tan- 
gential to the surface of the particle. Results are presented 
in Fig. [3j where the case of a passive nematic is shown as 
well for comparison, and where 77eff is plotted as a func- 
tion of in a double- logarithmic scale. In agreement with 
previous results and with Stokes' law j2T], we flnd that 
microrheology experiments performed in passive nematics 
(C = 0) lead to the measurements of quantitatively differ- 
ent viscosity, according to whether the particle is dragged 
along or orthogonally to the average director orientation 
in the bulk (with the "orthogonal" viscosity about a factor 
of 2 larger than the "parallel" one), in agreement with the- 
ory [34] and experiments [35]. In active nematics, contrac- 
tility increases the drag, while extensile activity reduces 
it - this is in line with the macrorheology results. 

More strikingly, our results show that as soon as / 0, 
a completely different behaviour is observed regarding the 
dependence of 77eff on R: Stokes' law is no longer valid 
when activity is turned on. Let us focus on Fig.|3](a) flrst: 
here 77eff rapidly increases with R when the particle is 
dragged along the director. When an extensile active fluid 
is considered (Fig. [3] (b)) instead, r^eff decreases with i?, 
the effect being more visible when the particle is dragged 
normal to the far fleld director. 

We qualitatively explain these results with the fol- 
lowing heuristic argument, sketched in Fig. [4] (and sum- 
marised previously in [16 ). If the particle were not mov- 
ing, the deformations in the director fleld due to the an- 
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Fig. 4. Sketch of the splay and bend deformation when the colloidal particle is dragged along different directions. Panels (a) 
and (b) refer to dragging in a contractile and extensile fluid respectively. Blue arrows refer to the active forces, while the red 
one represents the external force. The sketch is further discussed in the text. 
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choring condition would be perfectly symmetric, and forces 
on the probe due to the active stress would cancel. On the 
other hand, we expect that, when the particle is dragged 
through the active nematic fluid some asymmetries will 
show up in the director profile. In particular, in the case 
of planar anchoring, when the particle is dragged along 
the director, the splay in front of the particle should be 
greater than the one at the back (see Fig.|4^). This gener- 
ates a net force that opposes the external one, if the fluid is 
contractile, and that favours it, in the extensile case. That 
the effect is more important in contractile fluids can prob- 
ably be ascribed to the fact that they tend to splay more 
easily than their passive counterparts, while extensile ones 
are more stable in resisting splay [36]. An R dependence 
may then be expected as the splayed region depends on 
probe size. On the other hand, the director field bends, 
both down and upstream of the particle, when the exter- 
nal force is applied perpendicular to the far field director 
(see Fig. [4]3). When pulling upwards along that direction, 
again one expects the bending deformation on the top to 
be larger than the one at the bottom of the particle. These 
contributions combine to give a force opposing motion in 
contractile fluids and favouring it in extensile ones. Simi- 
larly to the case sketched in Fig. [4^, the effect is now larger 
for extensile fluids, as they tend to bend more easily than 
their passive counterpart, while contractile nematics are 
more resistant to bending [36]. 

This argument explains qualitatively our results on 
planar anchoring in Fig.[3| It further suggests that the ori- 
entation of the director field at the particle surface should 
play a key role in determining the microrheological drag 
of the particles. We therefore now investigate this aspect 
more in depth. We present in what follows both the cases 
of normal and of no imposed anchoring {W = 0). Although 
the latter case may be unlikely realised in experiments, 
where the colloidal probe is unlikely to not affect the ori- 
entation of the active nematogens nearby, this calculation 
leads to a useful limiting case. 

Results for normal anchoring are presented in Fig. [5j 
While pulling orthogonally to the far field director (see 
Fig. ^) leads to results which are broadly similar to the 
tangential anchoring case previously shown in Fig. [3) drag- 
ging along the director leads to a very different scenario. 
Now, quite strikingly, the particle moves much faster in 
contractile than it does in extensile nematics, as is ap- 
parent from Fig. ^jp. Remarkably, this can once more be 
explained by analysing the director deformation due to the 
anchoring condition, which is sketched in Fig. [6| In this 
case the splay is reversed with respect to the case of planar 
anchoring and one would expect to see an opposite effect, 
with contractile stresses pulling the particle forward, as 
we indeed observe. 

The scenario suggested by the cartoon in Fig.[6|is fully 
confirmed by a detailed analysis of the flow profiles (shown 
in Fig.[7|), around a probe particle being dragged through 
either passive nematic, a contractile active fluid, or an 
extensile one, and either along or perpendicular to the 
far field director. In particular, comparing panels (b), (d) 
and (f) one may easily note that in the contractile case the 



active stress leads to a flow field upstream of the particle 
which pulls it forward, hence facilitating its motion. In 
the extensile case the flow is smaller than in the passive 
case, again in agreement with our argument based on the 
director deformation. 

The result for normal anchoring, that the resistance to 
motion may now decrease in a contractile active fluid, is 
particularly important. Indeed it clearly shows that not 
only does the viscosity inferred from microrheological ex- 
periments strongly depend on the anchoring condition, 
but also that the very decrease/increase of the resistance 
to motion is not an inherent property of the underlying 
fluid. This is because the estimate of 77eff is also affected in 
a nontrivial way by the anchoring conditions, which will 
depend ultimately on the microscopic details of the inter- 
action between the colloidal probe and the particles that 
make up the active fluid. 

Intriguingly, unlike the other cases considered thus far, 
when pulling along the director in an extensile active fluid 
we find that there is a limiting force which needs to be 
overcome before the particle moves in steady state (see 
the inset in Fig. [Sf), where all curves refer to forces below 
this limiting value, so that ^ in steady state). This 
'yield force' may be seen as a microrheological analog of 
the yield stress in bulk rheology experiments. Its presence 
is equivalent to the absence of a linear regime which is 
why we do not report estimates for T^eff for this pulling 
mode. 

Finally, we turn to the case of no anchoring, W = 0. In 
this limit the director field is not affected by the presence 
of the particle, at least when the latter is quiescent. (When 
the particle moves the flow field will of course deform the 
nearby orientation of the active fluid.) This case might be 
thought of as more directly probing the actual properties 
of the fluid, as anchoring effects play no role here, and 
there is only a hydrodynamic coupling between probe and 
fluid. 

The resulting effective viscosity curves are shown in 
Fig. |8] They are qualitatively similar to those observed 
for planar anchoring (see Fig.|3|: contractility determines 
an increase in the measured viscosity with respect to the 
passive case, while a decrease is observed in extensile ac- 
tive nematics. However, no significant R dependence is 
found, presumably because now the director deformation 
is much smaller than in the W ^ case and is less affected 
by the probe size. 

All the results on microrheology presented so far were 
obtained with periodic boundary conditions. One may won- 
der whether this choice may lead to artifacts. We therefore 
also consider the case in which the system is bounded by 
two parallel planar walls. (This provides a physical way 
to stabilise the quiescent phase at low activity [Tl[2j[3l[4] 
although periodic boundary conditions also do this.) 

The results in this case are presented in Fig. |9] [37] . 
Pleasingly, it can be seen that they confirm, at least qual- 
itatively, what was previously found with periodic bound- 
ary conditions. The effect of activity is now slightly smaller, 
probably because the system is half the size of the one sim- 
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Fig. 5. (a) Dependence of the measured effective viscosity on radius for normal anchoring of the director field on the colloidal 
probe surface, at various C (see legend). The dragging direction is along the far field director. Other parameters are as in Fig.js] 
(b) Velocity versus time for F = 0.12, = 9.3 in the case of dragging along the far field director. The symbols are chosen 
according to the legend shown in (a). The inset in (b) shows the velocity of a particle of radius R = 9.3 pulled through an 
extensile gel (C = 0.0001) as a function of time for three different values of the extenal force F. Note the presence of a non-zero 
'yield force' further discussed in the text (all the curves shown in the inset correspond to forces before the yield force, so that 
V ^ ai late times). As in the previous cases, we used (periodic) cubic simulation box of volume V = 128^. 




Fig. 6. Sketch of the director field deformation for normal anchoring, and for a probe pulled along the far field director in a 
contractile fluid. As before, blue arrows refer to the active forces, while the red one represents the external force. 



ulated for the data shown in Fig. |3| so that the director 
deformation close to the particles may be slightly less. 

One important constraint which our model should ful- 
fill to reliably address an experimentally realisable situa- 
tion is that the probe used to obtain the microrheology 
data should be considerably larger than active nemato- 
gens. This is likely to be the case in practice if actomyosin 
(or some other molecular active fluid) is chosen as a host 
medium. If bacterial suspensions are considered, this con- 
dition would not easily be met, as the size of probes used 
in microrheology experiments is usually between 1 and 
10/im [12j, therefore comparable to the one of bacteria 
{E.coli for example is about ~l/im in size). If a sufficiently 
dense suspension is considered, the continuous approxima- 
tion still holds and our model still gives a good description 
of the system. However, one may imagine that in this case, 
rather than imposing any particular anchoring for the bac- 
terial director on its surface, the probe particle presence 
may simply enhance disorder in the bacterial orientation 
locally. One way to model this is to impose a null value of 
the amplitude of order 5'o for Q in eq. |3] The resultant di- 
rector field in the proximity of the particle is very similar 



to the W = case, as expected. The microrheology re- 
sults are shown in Fig. [lO] r^eff is again found to be larger 
with respect to the passive one in contractile nematics, 
and smaller in extensile ones, in agreement with most of 
our previous results. A rather mild dependence on R is 
observed for pulling both parallel and orthogonally to the 
far-field director and, remarkably, no significant difference 
is observed between these two cases in terms of the R de- 
pendence with this boundary condition. This is probably 
because the chosen boundary condition melts the active 
nematic close to the particle surface, so that the defor- 
mations of the orientational order close to the particles 
should not much depend on the dragging direction. 



4 Discussion and conclusion 

In conclusion, we presented numerical results on the bulk 
rheology and on the microrheology of active fluids, specif- 
ically active nematics, consistent with previous work [2j 
[3l[4l[9]. We showed that when shearing an active nematic 
fluid in a quasi-one dimensional geometry, the effective 
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Fig. 7. Flow field patterns close to a moving particle of radius R— 11.3 subject to a force F — 0.12, in a passive nematic liquid 
crystal (a,b), in a contractile active fluid (c,d), and in an extensile active fluid (e,f). In all cases there is normal anchoring of the 
director field at the particle surface. In (a,c,e) the pulling direction is perpendicular to the far field director, whereas in (b,d,f) 
the colloid is dragged par ailed to the far field. 
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Fig. 8. Dependence of the effective viscosity on particle size for no anchoring, and different C (see legend), (a) and (b) describe 
a contractile and extensile active fluid respectively. As previously, filled and open symbols correspond to pulling along and 
perpendicular to the far field director respectively and we used (periodic) cubic simulation box of volume V = 128"^. 
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homogeneous anchoring of the director field at the wall of the container. The anchoring determines in this case the far field 
director profile. All parameters are as in Fig. [3] apart from the system size which is 64 along the direction orthogonal to the 
walls, and 128 in the other two directions. 



2.7 



viscosity strongly depends on system size, and also on the 
imposed anchoring of the active nematogen orientation at 
the wall. Our numerical results are in very good agreement 
with analytical estimates, obtained by linearising the rele- 
vant hydrodynamic equations of motion. Our results thus 
demonstrate that bulk rheology experiments performed 
on active nematics should strongly depend on the details 
of the measurement setup, precluding any identification 
of the effective viscosity of an active fluid as an inherent 
property of the material. 

We also simulated the simplest possible microrheol- 
ogy experiment in an active fluid. That is, we applied 
a constant force, either parallel or orthogonal to the far 
field director, to a spherical probe particle embedded in 
an active nematic, and measured the steady state veloc- 
ity it attains, which leads to an estimate for the drag on 
the particle. If Stokes' law held in this active context, we 
would expect the drag coefficient to be linearly propor- 



tional to the particle size and to an effective viscosity, 
77eff- We presented results for both planar and normal an- 
choring, which show that the colloidal drag is in both cases 
strongly non-Stokesian, as the dependence on size is highly 
non-linear. We gave a qualitative explanation of this non- 
Stokesian drag in terms of the local deformation of the 
director orientation close to the particle surface. Our ar- 
gument also explains why planar and normal anchoring 
lead to qualitatively different results: for example, when 
a colloidal probe is pulled along the far field director, its 
drag is larger than that in a passive nematic when the an- 
choring is planar, but it decreases for normal anchoring. 

We also presented numerical results for the effect on 
microrheology experiments, with different boundary con- 
ditions. In particular, we simulated the case in which the 
sample is bounded by solid walls, as opposed to that in 
which periodic boundary conditions are considered. Pleas- 
ingly, we found that our results in the two cases are sim- 



G. Foffano et al.: Bulk rheology and microrheology of active fluids 



11 




_ f 


= 




' 1 ' 




= 0.0001 






» — 


e 9- — ®"~o--a> 














a-^ 


















1 




1 


1 1 1 



1.5 



InR 



2.5 



Fig. 10. Dependence of the effective viscosity on particle size with Q = on its surface, and different values of C (see legend), 
(a) and (b) refer to a contractile and extensile active fluid respectively. As previously, filled and open symbols correspond to 
pulling along and perpendicular to the far field director respectively and we used (periodic) cubic simulation box of volume 
V = 128^ 



ilar. Finally we also considered both an idealised case in 
which the colloidal probe does not affect the orientation 
of the active nematogens in equilibrium, and the extreme 
case in which instead the order parameter is fixed to zero 
on its surface. With both boundary conditions, we found 
that the results qualitatively resembled those obtained 
with planar anchoring, and we observed that the drag on 
the particle, hence r^eff, increases in contractile fluids, and 
decreases in extensile ones, although the dependence on 
size was milder. 

We hope our results will stimulate further experimen- 
tal work, on both the bulk rheology and the microrheology 
of active fluids such as actomyosin solutions or very dense 
bacterial suspensions. (The latter can be probed via opti- 
cal tweezers [38 .) 
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